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Abstract 

A novel approach to parton distributions parameterization in terms of quantum 
statistical functions is here outlined. The description, already proposed in pre- 
vious publications, is here improved by adding to the statistical distributions an 
unpolarized liquid component. This new contribution to fermion partons is able to 
reproduce the expected low x behaviour of structure functions. The analysis pro- 
vides a satisfactory description of polarized and unpolarized deep inelastic data and 
shows a possible connection between the Gottfried and Bjorken sum rules. 
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1 Introduction 



In a recent series of papers the role of Pauli exclusion principle to explain the 

experimental data on unpolarized and polarized structure functions of the nucleons has 
been studied. 

The relevance of the quantum statistics on the parton distribution is supported by 
several phenomenological observations. The most relevant phenomenon is certainly the 
measurement of a defect in the Gottfried sum rule [l^. It can be explained in terms 
of a Pauli blocking effect on the production of u sea quark with respect to d sea quark, 
which yields a flavour asymmetry between u and d in the proton. Another interesting 
observation is a relationship which seems to occur between the shapes of the quark parton 
distributions and their first momenta, which is the typical characteristic of Fermi-Dirac 
distribution functions. 

These considerations naturally suggest to use quantum statistically inspired param- 
eterizations for the parton distribution functions. In this description, the independent 
variable which plays the role of the energy is the Bjorken variable x, and the distribu- 
tions assumed will be of the Fermi-Dirac kind for quarks, and Bose-Einstein for gluons. 
Interestingly, to properly describe the low x behaviour of the structure functions, a liquid 
unpolarized component dominating the very low x region has to be added. It does not 
affect the quark parton model sum rules (QPMSR) but it is necessary to reproduce the 
antiquark distribution at low x. 

The paper is organized as follow, in Section 2 we summarize the phenomenological 
motivations behind our description for the parton distributions. In Section 3 the distri- 
bution function parameterizations are shown in detail and discussed in connection with 
a possible physical interpretation. Section 4 deals with the results of a fitting procedure 
performed to get the free parameters of the distribution functions. The theoretical pre- 
dictions are shown in comparison with the experimental data and the results for QPMSR 
are also discussed. In Section 5 we give our conclusions and remarks. 
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2 Evidence for quantum statistical effects in parton 
distributions 



As already stated, an experimental evidence for a central role played by the Pauli principle 
in the physics of nucleon is the defect in the Gottfried sum rule 

"1 1 . ™. . 1 



I -M{x) - F^{x)] dx = -{u + u-d~d) , (1) 
Jo X 3 



that for SU{2)i invariant sea quark distributions {d = u) gives Iq = 1/3. Indeed NMC 
experiment |TD[ measures for the l.h.s. of Eq. (|T|) 

Ig = 0.235 ± 0.026 , (2) 

implying 

J- M = 0.15 ±0.04 u-d = 0.85 ±0M . (3) 

This experimental result can be explained following the idea of Field and Feynman [p 



who suggested that, in the proton, the Pauli principle depresses the production of uu 
pairs in the proton with respect to dd, since it contains two valence u quarks and only one 
d. We will return on this fact in the following to focus on a possible connection between 
the violation of the above sum rule and the Bjorken one |]T2| . 



From the previous considerations one expects a relevant role played by the statistics 
in the whole phenomenology of deep inelastic scattering, and thus it suggests to look for 
others typical characteristics of this behaviour. 

A peculiar characteristic of a Fermi-Dirac statistical function is certainly the strong 
connection between shape and abundance. This is an immediate consequence of the Pauli 
exclusion principle which forbids two or more fermions to have the same quantum numbers 
and implies that the more abundant is distribution the broader is in x the associated 
function. 

With the aim to check if this situation occurs in the nucleon structure let us consider 
the abundances of valence quarks in the nucleons. As it is well-known, at = they 
are connected to the axial couplings of the baryon octet F and D, through the relations 

«L = i + ^ «L = i-^ , (4) 

d^ -l±l^ dl _ ^-F + D 



By using the experimental values obtained by the two bodies strong decays of hyperons 
F = 0.464 ± 0.009 D = 0.793 ± 0.009 H, fg, we get 



T i 



+ dlai + dL , (6) 



which tells us that ul^i is the most abundant parton in the proton, at least at = 0. 
Moreover, by observing that F ~ 1/2 and D ~ 3/4 one also gets 

1 _ 1 _ dial + dial /yN 
"'val ~ 2 ~ 2 ' 

In addition to this, the behaviour at high x of F2"(x)/Ff (x) known since a long time. 



and the more recent polarization experiments [|TB[, |T^, which show that at high x the 
partons have spin parallel to the one of the proton, imply that u\x) is the dominating 
parton distribution in the proton at high x. Thus, to the most abundant corresponds 
effectively a broader distribution in the Bjorken variable x. 

Eq. has also another interesting implication. In a previous work we assumed that 
the parton distributions at a given large depend on their first momenta (abundances) 
computed at = 

p{x) = T{x,pvai) , (8) 

with J-" an increasing function of p^ai and with a broader shape for higher values of p^ai- 
From this assumption and by virtue of we get 
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u^(x) = - [d^{x) + d^{x)\ = -d{x) , (9) 
which implies 

Au{x) = u\x) — (x) = u{x) — d{x) . (10) 

Note that, Eq. (|TUp connects the contribution of Au{x) to gi{x), with the terms due to 
up and down quarks in the unpolarized structure functions of nucleons F2{x) 

^9Ux)L = l[Fiix)-F-{x)l^, . (11) 

This relation should hold in good approximation for the total quantities xgi{x) and 
(a:) — -F^(x), since the contribution in (?i(x) due to Ad{x) is depressed for the twofold 
reason that = (l/4)e^ and Ad^ai — —{^/4:)Auvai- By integrating Eq. (|T^) one thus get 
a connection between the spin sum rule and the Gottfried sum rule and in turn a relation 
between their possible defects. 



3 Parton distributions as Fermi— Dirac and Bose— Einstein 
statistical functions 



The previous considerations on the role played by the Pauli principle in the nucleon 
structure suggest to assume Fermi-Dirac distributions in the variable x, at least for large 
X, for the quark partons 



Px{x) = f{x) 



( X - xi\)\ 



-1 



'12) 



where the index A denotes the different species of quarks, characterized by flavour and 
polarization. In Eq. ([T^), f{x) is a weight function which accounts for the energy level 
density, and because it is connected to the nonperturbative aspect of QCD results indepen- 
dent of flavours and polarization. The universal parameter x represents the temperature 
for the system, whereas x{\) stands the thermodynamical potential of the parton A. 

The expression chosen for f{x) is inspired by the expected power behaviour at x = 0, 
and by the obvious kinematical cut which forces the function to vanish at a; = 1. In order 
to satisfy these constraints we assume for simplicity a power low dependence on x 

f{x) = Ax'^il-xf . (13) 



Indeed, the assumption that the form given in Eq. (|1^) for the quark distribution func- 
tions, which requires different thermodynamical potentials in order to describe the ex- 
perimental data, is valid in the low x limit as well has at least two unpleasant features. 
Firstly, one gets in the nonperturbative region different behaviour for the different par- 
ton distributions, where on the contrary one would expect an universal dependence on 
X. Moreover, the power dependence on x of Eq. ([T^), fitted by the experimental data 
mostly placed at the large x, is not suitable to reproduce the more divergent contribution 
expected at low x. This most divergent part does not contribute to QPMSR as the ones 
given by Gottfried and Bjorken with 1 = 1 quantum numbers exchanged. To this aim 



we add to (|12[) an unpolarized component, which we call liquid to stress the possibility 
that it is connected to the presence, at low x, of a new phase in the quark-gluon plasma 
due to the highly nonperturbative QCD regime 



p^[x) = — x"^(l - xf^ + A x"(l - xf 



fx — x(A) \ 



-1 

(14) 



In the fitting procedure we take as free parameters, apart from the constants involved in 
f(x) and in the liquid component of (Q), the temperature x and the x for u^^^\ d^^-''\ u and 
d (the latters are assumed not polarized). We tried initially to introduce spin- dependent 
x's also for the g's and to test the relationship 

Au(x) = u{x) — d{x) , (15) 

assumed in previous works @, ||, but unfortunately, we found practically the same with 
negative and positive values for A'u(x) and/or Ad{x). Hence, for not loosing predictivity 
in the fit procedure we have assumed unpolarized antiquarks. 

As far as the strange quarks are concerned, we assume for simplicity unpolarized 
distribution functions given by the empirical expression 

. = ^-=^" . (16) 

which experimentally is very well satisfied. 

Analogously, for the gluons, if we neglect their polarization , the bosonic statistic 
suggests the consider, at large x, the distribution function 

16 r ^ X n-1 



G(x) = —A a;"(l - xf [exp ( 

3 I \ X / 



1 



(17) 



where the factor 16/3 is just the product of 2 (S'^(G) = ±1) times 8/3, the ratio of the 
colour degeneracies for gluons and quarks. 

4 Discussion of the results 

By assuming for the parton distributions Eqs. (p^ and (plT]), we fit the distribution 
parameters from the experimental data for -Ff (x) — F^(x) [|10|, xF^i^x) [1^, xg^x) [|1^, ^ 



and xgi{x) [|T9|, which do not receive contributions from the liquid component, and from 



F^{x)/Fi{x) [ig anda;g(a;) |2g . 

The available experimental data on deep inelastic scattering observables correspond 
in general to different values of Q^. This would suggest, in order to use the data to 
determine the distribution parameters, which in general will depend on Q^, to apply the 
evolution equations to lead all the experimental results to the same Q^. In our analysis 
we have neglected this dependence of the distribution parameters, since we expect 



from the evolution equations a smooth logarithmic dependence. As far as the polarized 
distributions is concerned, in fact, the expected dependance results to be smaller 
than the experimental errors, and thus our analysis is slightly affected by neglecting this 
dependance. 

Indeed, the data on unpolarized nucleons structure functions are at = 4 GeV^ , 
the neutrino data at = 3 GeV^ []rB[, and q measures are performed at = 3 GeV^ 
and 5 GeV^ |20] and differ at small x, while our curve is intermediate between the two 
sets of data. The data on gi{x) are at Q"^ = 2 GeV^ [|l9l, whereas gi{x) is measured at 
Q2 = 10 GeV^ by SMC [|6l and at = 3 CeV^ by E143 despite some narrowing 
of the distribution at higher showing up in the data, the values of Ip are in good 
agreement. 

In Table 1 we report the parameters found in the present analysis [Qand compare them 
with the results of a previous fit (without liquid) 0, and with the ones by Bourrely and 
Soffer found on similar principles, but with several different assumptions. In the Figures 
1.-6., the predictions for the nucleon structure -^^(a;) — ^^(a;)/^^^), xgi{x), 

xgi{x), xFslx) and xq{x) are shown, respectively, and compared with the experimental 
data. 

The parton distributions found in |Q are described in Figure 7. Since the total momen- 
tum carried by fermion partons is 53%, we get Xq = —1/15 by requiring that the gluons 
carry out the remaining part of the proton momentum. In Ref . , the gluon distribution 



is compared with the information found on them in CDHSW |21|], SLAC+BCDMS 
and in NMC experiments at = 20 GeV^. The agreement is fair for x > .1, while 



the fast increase at small x, confirmed also from the data at very small x at Hera [|4 



confirms that a liquid component is needed also for gluons. The excess at high x of our 
curve with respect to experiment may be, at least in part, explained by the expected nar- 
rowing of the distribution from = 4 GeV^, where we fit the unpolarized distributions, 
to Q2 = 20 GeV\ 

The inclusion of the liquid term and the extension of our fit to the precise experimental 
results on neutrinos has brought to substantial changes in the parameters with respect 
to the previous work 

The low X behaviour of f{x) become smoother (~ 3;--203±.oi3 instead of x~^'^^), but this 
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is easily understood since the previous behaviour was a compromise between the smooth 
gas component and the rapidly changing liquid one to reproduce the behaviour of q{x). 
The liquid component, relevant only at small x, carries only .6% of parton momentum 
and its behaviour ~ x~^'^^, similar to the result found in |2^, is less singular than the 
one, suggested in the framework of the multipherial approach to deep-inelastic scattering, 
proportional to ~ x~^'^ [EBf. The parameter x{u^) took the highest value allowed by us 



(1.), since the factor in /(x), (1 — x)^'^^, is taking care to decrease u\x) at high x. The 
temperature x is larger than the previous one and the one found by Bourrely and Soffer 
0. Instead x{u^) is slightly smaller than the previous determination |^ and about half 
the value found in 0, where f{x) is different for u'^ and uK 

The ratio r = u^{x)/d{x) varies in the narrow range (.546, .564) in fair agreement with 
the constant value 1 — F = .536 ± .009 assumed in and slightly larger than the value 
1/2 taken in and §. 

The central value found for the first moment of Ugas{x), -03, is smaller than dgas{x)/2, 
.08, while Eq. (p!5| ) implies u{x) > d{x)/2. However, the large upper error on Ugas and 
the uncertainty in disentangling the gas and liquid contributions for the g's do not allow 
to reach a definite conclusion about the validity of Eq. (|T5|). 

Finally, we can compare the predictions of with the measured asymmetry for Drell- 
Yan production of muons at y = in pp and pn reactions 

^ d(Tpp/ dy - daprj dy ^^^^ 
dapp/dy + dapn/dy 

which at rapidity y = reads 

_ jX^jx) - l)(4A(x) - 1) + iXjx) - l)(4A,(x) - 1) 
"""^ (A,(x) + l)(4A(x) + l) + (A(a;) + l)(4A,(x) + l) ' ^ ^ 

where Xs{x) = u{x)/d{x) and A(x) = u{x)/d{x). At x = .18 we have As(.18) = .454 and 
A (.18) = 1.748 giving rise to A£,y{.18) = —.138 in fair agreement with the experimental 
result -.09 ±.02 ±.025 |3 . 



The behaviour of A£)y{x) is plotted in Figure 8 together with the experimental point 
measured by NA51 collaboration. 
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5 Conclusions 



We compared with data the quark-parton distributions given by the sum of Fermi-Dirac 
functions and of a term not contributing to the QPMSR relevant at small x. We obtain 
a fair description for the unpolarized and polarized structure functions of the nucleons as 
well as for the F^lx) precisely measured in (anti)neutrino induced deep- inelastic reactions 
and for q total distribution. The conjectures of previous works on d distributions are 
well confirmed by the values chosen for their thermodynamical potentials. As long as the 
implications for QPMSR the values found for the first momenta of the various parton 
species give l.h.s.'s consistent with experiment. For the fundamental issue of the Bjorken 
sum rule, as advocated in previous works |jl|, ^ and f^, we get 

Au ~ u - d + 2F -1 , (20) 
Ad > F-D , (21) 

to confirm the suspicion of a violation of Bjorken sum rule related to the defect in the 
Gottfried sum rule. 

A word of caution is welcome for our conclusions on the violation of Bj sum rule, since 
we did not include the effect of QCD corrections in relating the quark parton distributions 
to the structure functions. Also we assumed no polarization for q, being unable to get a 
reliable evaluation of Aq with the present precision for the polarized structure functions 
at small x. Indeed our description of gf{x) and gi{x) is good in terms of Au{x) and 
Ad{x), but our prediction is smaller than the central values of the three lowest x values 
measured by SMC. 
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Parameters 


Previous fit ||5| 


Fit BS § 


Actual fit P 










= 2.47 


A 


.58 




+ 09 


a 


-.85 


-.646 for ul^i 
-.262 for , 


-.203 ± .013 








o +-05 
2-34 _.06 


A 

Al 






0895 +-°^°^ 
■^^^^ -.0084 


OIL 






-1.19 ± .02 


Q 

PL 






^•66 -1.59 


X 


.132 


.092 


.235 ± .009 


gas abund. 


x{v)) 


.524 


.510 for ul^i 


1.00 ±.07 


1.15 ± .01 


x{u^) 


.143 


.231 for uL, 


.123 ± .012 


.53± .01 








-.068 _ 


.33± .03 


x{S) 






200 


.62± .01 


x[u^) 


—.216 




-.886 ± .266 


015 

•^^^ -.009 


x{u^) 


-.141 




// 


// 


x{d}) = 


// 




- 460 +-°^' 
-.064 


08 
.08 



Table 1. Comparison of the values for the parameters of our best fit |^ with the corresponding 
quantities, if any, found in previous analysis (|], Q- 



Sum rule 



Experimental data 



Our fit 



QPM 



GLS 
G 



EJ 



Bj 



2. 50 ±.018 ±.078 [ig 



.235 ±.026 10 



.136 ±.011 ±.011 [|T§ 
.129 ± .004 ± .009 H 

-.022 ± .011 m 



2.44 



-.04 



.07 
.20 ±.02 

.122 ± .007 
-.030 ±.010 

.152 ± .010 



3 

1/3 
.188 ± .005 

-.021 ± .005 
.209 



Table 2. Comparison of our predictions for the sum rules with the experimental 
values and with the quark parton model (QPM) predictions without QCD corrections. 
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Figure 1: The prediction for (x) — is plotted and compared with the experimental 

data [10] 
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Figure 4: xgi{x) is plotted and compared witli tlie data [21]. 
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